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Abstract 

We explore in the framework of Quantum Computation the notion of Computabil- 
ity, which holds a central position in Mathematics and Theoretical Computer Sci- 
ence. A quantum algorithm for Hilbert's tenth problem, which is equivalent to the 
Turing halting problem and is known to be mathematically noncomputable, is pro- 
posed where quantum continuous variables and quantum adiabatic evolution are 
employed. If this algorithm could be physically implemented, as much as it is valid 
in principle — that is, if certain hamiltonian and its ground state can be physically 
constructed according to the proposal — quantum computability would surpass clas- 
sical computability as delimited by the Church- Turing thesis. It is thus argued that 
computability, and with it the limits of Mathematics, ought to be determined not 
solely by Mathematics itself but also by Physical Principles. 



1 Introduction 

Computation based on the principles of Quantum Mechanics [1] has been shown to offer 
better performances over classical computation, ranging from the square-root improvement 
in an unstructured search [2] to the exponential gain in the factorisation of integers [3]. 
However superior in reducing the complexity of hard computation, these quantum algo- 
rithms and all the others discovered so far are only applicable to the classically computable 
functions. That leaves untouched the class of classically noncomputable functions, such as 
the halting problem for Turing machines [4]. It is in fact widely believed that quantum 
computation cannot offer anything new about computability [5]. Contrary to this, we pro- 
pose that quantum computation may be able to compute the noncomputables, provided 
certain hamiltonian and its ground state can be physically constructed. We propose a 
quantum algorithm for the classically noncomputable Hilbert's tenth problem [6] which 
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ultimately links to the halting problem for Turing machines in the computation of partial 
recursive functions. 

The practical details of implementation the quantum algorithm for this class of problems 
are not considered in this conceptual study, and will be investigated elsewhere. 

2 Hilbert's tenth problem 

At the turn of the last century, David Hilbert listed 23 important problems, among which 
the problem number ten could be rephrased as: 

Given any polynomial equation with any number of unknowns and with in- 
teger coefficients: To devise a universal process according to which it can be 
determined by a finite number of operations whether the equation has integer 
solutions. 

This decision problem for such polynomial equations, which are also known as Diophan- 
tine equations, has eventually been shown in 1970 by Matiyasevich to be undecidable [6, 7] 
in the Turing sense. It is consequently noncomputable/undecidable in the most general 
sense if one accepts, as almost everyone does, the Church- Turing thesis of computability. 
Since exponential Diophantine, with the unknowns in the exponents as in the example of 
the Fermat's last theorem, can be shown to be Diophantine with supplementary equations, 
the study of Diophantine equations essentially covers the class of partial recursive func- 
tions, which is at the foundations of classical algorithms. The undecidability result is thus 
singularly important: Hilbert's tenth problem could be solved if and only if the Turing 
halting problem could be. 

3 Turing halting problem 

The halting problem for Turing machines is also a manifestation of undecidability: a Turing 
computation is equivalent to the computation of a partial recursive function, which is only 
defined for a subset of the integers; as this domain is classically undecidable, one cannot 
always tell in advance whether the Turing machine will halt (that is, whether the input 
is in the domain of the partial recursive function) or not (when the input is not in the 
domain) . 

A version of the proof of the unsolvability of the halting problem based on the Cantor 
diagonal argument goes as follows. The proof is by contradiction with the assumption of 
the existence of a computable halting function h(p, i) which has two integer arguments - p 
is the Godel encoded integer number for the algorithm and i is its (encoded) integer input: 




(1) 
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One can then construct a program r(n) having one integer argument n in such a way that 
it calls the function h(n, n) and 

r{n) halts if h(n,n) = 1 

r(n) loops infinitely (i.e., never stops) otherwise. 
The application of the halting function h on the program r and input n results in 

h(rn) - { ° if/i (n,n) = l 

n[r ' n) ~ \ 1 if h(n,n) = [Z) 

A contradiction is clearly manifest once we put n = r in the last equation above. 

The construction of such program r is transparently possible, unless the existence of a 
computable h is wrongly assumed. Thus the contradiction discounts the assumption that 
there is a classically algorithmic way to determine whether any arbitrarily given program 
with arbitrary input will halt or not. 

This contradiction argument might be side stepped if we distinguish and separate the 
two classes of quantum and classical algorithms. A quantum function qh(p,i), similar to 
eq. (1), can conceivably exist to determine whether any classical program p will halt on 
any classical input % or not. The contradiction in eq. (2) would be avoided if the quantum 
halting qh cannot take as argument the modified program r, which is now of quantum 
character because it now has quantum qh as a subroutine. This will be the case if qh can 
only accept integer while quantum algorithms, with proper definitions, cannot in general 
be themselves encoded as integers. It is clear even in the case of a single qubit that the 
state a|0) cannot be encoded as integers for all a and /3 - simply because of different 

cardinalities. In fact, the no-cloning theorem [8] of quantum mechanics does restrict the 
type of operations available to quantum algorithms. 

In essence, the way we will break the self-referential reasoning here by the differentiation 
between quantum and classical algorithms is similar to the way John von Neumann and 
Bertrand Russell resolved the set theory paradox (to do with "The set of all sets which 
are not members of themselves" ) by the introduction of classes as distinct from sets. (For 
other lines of arguments, see [9].) 

4 An observation 

It suffices to consider nonnegative solutions, if any, of a Diophantine equation. Let us 
consider the example 

(x+lf + (y + lf-(z + lf + cxyz = 0, ceZ, (3) 

with unknowns x, y, and z. To find out whether this equation has any nonnegative integer 
solution by quantum algorithms, it requires the realisation of a Fock space built out of the 
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"vacuum" |O a ) by repeating applications of the creation operators a x , a y and a\, similarly 
to that of the 3-D simple harmonic oscillators. 

[aj,a]] = l for j = x, y, z, (4) 
[a k , dj] = [a k , a]] = for j ^ k. 

Upon this Hilbert space, we construct the hamiltonian corresponding to (3) 

Hp = {{a x a x + l) 3 + (a y a y + l) 3 - (a\a z + l) 3 + c(a\a x )(a y a y )(a\a z )^ , 

which has a spectrum bounded from below - semidefinite, in fact. 

Note that the operators Nj = a jo,- have only nonnegative integer eigenvalues rij, and 
that [Nj, Hp] = = [N, Nj] so these observables are compatible - they are simultaneously 
measurable. The ground state \g) of the hamiltonian so constructed has the properties 

N i\9) = n i\9), 

H P \g) = ((n x + l) 3 + (n y + l) 3 -(n z + l) 3 + cn x n y n z ) 2 \g)=E g \g), 
for some (n x ,n y ,n z ). 

Thus a projective measurement of the energy E g of the ground state \g) will yield 
the answer for the decision problem: The Diophantine equation has at least one integer 
solution if and only if E g — 0, and has not otherwise. (If c = in our example, we know 
that Eg > from Fermat's last theorem.) 

If there is one unique solution then the projective measurements of the observables 
corresponding to the operators Nj will reveal the values of various unknowns. If there are 
many solutions, finitely or infinitely as in the case of x 2 + y 2 — z 2 = 0, the ground state \g) 
will be a linear superposition of states of the form \n x ) <8> \n y ) <E> \n z ), where (n x , n y , n z ) are 
the solutions. In such situation, the measurement may not yield all the solutions. However, 
finding all the solutions is not the aim of a decision procedure for this kind of problem. 

Notwithstanding this, measurements of Nj of the ground state would always yield some 
values (n x ,n y ,n z ) and a straightforward substitution would confirm if the equation has a 
solution or not. Thus the measurement on the ground state either of the energy, provided 
the zero point can be calibrated, or of the number operators will be sufficient to give the 
result for the decision problem. 

The quantum algorithm with the ground-state oracle is thus clear: 

1. Given a Diophantine equation with K unknowns x's 

D{ Xl ,---,x K ) = 0, (5) 
we need to simulate on some appropriate Fock space the quantum hamiltonian 

Hp = (p(a\a u --- ,a ] K a K )f . (6) 
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2. If the ground state could be obtained with high probability, measurements of appro- 
priate observables would provide the answer for our decision problem. 

The key ingredients are the availability of a countably infinite number of Fock states, 
the ability to construct /simulate a suitable hamiltonian and to obtain/identify its ground 
state via quantum measurements. As a counterpart of the semi-infinite tape of a Turing 
machine, the Fock space is employed here instead of the qubits of the more well-known 
model of quantum computation. Its advantage over the infinitely many qubits which would 
otherwise be required is obvious. 

5 Some preliminary comments 

We do not look for the zeroes of the polynomial, D(xi, • • • , which may not exist, but 
instead search for the absolute minimum of its square which exists, 

< min (D(x u • • • , x K )f < (D(0, • • • , 0)) 2 , 

and is finite because lim^oo (D(xi, • • • , xk)) 2 diverges. 

While it is equally hard to find either the zeroes or the absolute minimum in classical 
computation, we have converted the problem to the realisation of the ground state of 
a quantum hamiltonian and there is no known quantum principle against such act. In 
fact, there is no known physical principles against it. Let us consider the three laws 
of thermodynamics concerning energy conservation, entropy of closed systems and the 
unattainability of absolute zero temperature. The energy involved in our algorithm is finite, 
being the ground state energy of some hamiltonian. The entropy increase which ultimately 
connects to decoherence effects is a technical problem for all quantum computation in 
general, and we will discuss this further below. As we will never obtain the absolute 
zero temperature, we only need to satisfy ourselves that the required ground state can be 
achieved with a more-than-even chance. Then there is a probability boosting technique, 
see later, to bring that chance to as closed to unity as one pleases. 

It may appear that even the quantum process can only explore a finite domain in a 
finite time and is thus no better than a classical machine in terms of computability. But 
there is a crucial difference. 

In a classical search even if the global minimum is come across, it cannot generally be 
proved that it is the global minimum (unless it is a zero of the Diophantine equation). 
Armed only with mathematical logic, we would still have to compare it with all other 
numbers from the infinite domain yet to come, but we obviously can never complete this 
comparison in finite time -thus, mathematical noncomputability. 

In the quantum case, the global minimum is encoded in the ground state. Then, by 
energy tagging, the global minimum can be found in finite time and confirmed, if it is the 
ground state that is obtained at the end of the computation. And the ground state may 
be identified and/or verified by physical principles. These principles are over and above 
the mathematics which govern the logic of a classical machine and help differentiating 
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the quantum from the classical. Quantum mechanics could "explore" an infinite domain, 
but only in the sense that it can select, among an infinite number of states, one single 
state (or a subspace in case of degeneracy) to be identified as the ground state of some 
given hamiltonian (which is bounded from below). This "sorting" can be done because 
of energetic reason, which is a physical principle and is not available to mathematical 
computability. 

On the other hand, our proposal is apparently in contrast to the claim in [5] that quan- 
tum Turing machines compute exactly the same class of functions as do Turing machines, 
albeit perhaps more efficiently We could only offer here some speculations about this 
apparent discrepancy. The quantum Turing machine approach is a direct generalisation of 
that of the classical Turing machines but with qubits and some universal set of one-qubit 
and two-qubit unitary gates to build up, step by step, dimensionally larger, but still di- 
mensionally finite unitary operations. This universal set is chosen on its ability to evaluate 
any desirable classical logic function. Our approach, on the other hand, is from the start 
based on infinite-dimension hamiltonians acting on some Fock space and also based on the 
special properties and unique status of their ground states. The unitary operations are 
then followed as the Schrodinger time evolutions. Even at the hamiltonian level higher 
orders of the operators a and a', i.e. not just two-body but many-body interactions in a 
sense, are already present. This proliferation, which is even more pronounced at the level 
of the time-evolution operators, together with the infinite dimensionality and the unique 
energetic status of the vacuum could be the reasons behind the ability to compute, in a 
finite number of steps, what the dimensionally finite unitary operators of the standard 
quantum Turing computation cannot do in a finite number of steps. Note that it was the 
general hamiltonian computation that was discussed by Benioff and Feynman [10, 11] in 
the conception days of quantum computation. 

Indeed, Nielsen [12] has also found no logical contradiction in applying the most general 
quantum mechanical principles to the computation of the classical noncomputable, unless 
certain Hermitean operators cannot somehow be realised as observables or certain unitary 
processes cannot somehow be admitted as quantum dynamics. And up to now we do 
not have any evidence nor any principles that prohibit these kinds of observables and 
dynamics. (Ozawa [13] has produced some counter arguments but we think they are not 
quite applicable here. See [9].) 

Our general algorithm above could be realised by, but in no way restricted to, the 
following methods to simulate the required hamiltonian and to obtain the ground state 
adiabatically. 

6 Simulating the hamiltonians 

One way to construct any suitable hamiltonian so desired is through the technique of 
ref. [14]. We consider the hermitean operators, where j is the index of the unknowns of 
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the Diophantine equation, 



Xj = -^(cij + aj), 

Pj = -^(aj-at), (7) 

Together with the availability of the fundamental hamiltonians 

X v Pj, (X 2 + P 2 ), ±(X k P j + P,X fc ), and (X* + P 2 ) 2 (8) 

one could construct the unitary time evolutions corresponding to hamiltonians of arbitrary 
hermitean polynomials in {Xj, Pj}, and hence in {a]aj}, to an arbitrary degree of accuracy. 
These fundamental hamiltonians correspond to translations, phase shifts, squeezers, beam 
splitters and Kerr nonlinearity. 

With the polynomial hamiltonian constructed, we need to obtain its ground state. Any 
approach that allows us to access the ground state will suffice. One way is perhaps to 
use that of quantum annealing or cooling [15]. Another way is to employ the quantum 
computation method of adiabatic evolution [16]. 



7 Adiabatic quantum evolution 

In the adiabatic approach, one starts with a hamiltonian Hj whose ground state \gj) is 
readily achievable. Then one forms the slowly varying hamiltonian 7i(s), s — -f, which 
interpolates between Hj and H P in the time interval t G [0, T] 

H(s) = (l-s)H I + sH P . (9) 

Note that we can replace this linear interpolation by some non-linear one provided the 
conditions of the adiabatic theorem are observed. According to this theorem, the initial 
ground state will evolve into our desirable ground state \g) up to a phase: 

Urn Texp{-iT f 1 n(r)dr}\ gi ) = e*\g). (10) 

T^oo I JO J 

For the hamiltonian (9), an estimate of the time T after which the system remains with 
high probability in the ground state is 

T » 11 Hl ~ 2 Hp K (11) 



with 



Hj-HpW = max | (e(t) | (H T - H P ) \g(t))\ , (12) 
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and 



mm (E e (t) - E g (t)) , (13) 



where \g(t)) and \e(t)) are respectively the instantaneous ground state and the first excited 
state of (9) with instantaneous eigenvalues E g (t), E e (t). 

The time-ordering operator on the left hand side of (10) can be approximated as 

exp{-iTH(T N )Ar} ■ ■ • exp{-iTH(r 1 ) At} 

for [16] 

NAt = 1, 

At II Hi - Hp || < 1. (14) 

Note that we have employed here the "norm" || . || as defined in (12) for the various 
hamiltonians which are unbounded from above. This norm is the relevant measure for 
the problem is only concerned with the lowest states of the interpolating hamiltonian (9). 
In each interval Ar, the unitary operators exp{— iH(r k )TAT}, for k — 1, • • • , JV, can be 
expressed through the subdivision of TAr into m subintervals of sufficiently small size 5s 
satisfying m5s = TAt, 

exp{-iH(r k )TAr} = (exp{-iH(T k )5s}) m . 

Each of the m factors on the right hand side of the last expression can be now simulated 
through the approach of [14], where it was shown that the number of steps M grows as a 
small polynomial in the order of the polynomial in the hamiltonian 7i{r k ) to be simulated, 
the accuracy to be enacted, and the time interval TAr over which it is to be applied. 

In this way, the requirements of the adiabatic conditions on the one hand and of, on the 
other hand, the simulations of the hamiltonians in the time interval TAr can be satisfied. 



8 An adiabatic algorithm 

In order to solve Hilbert's tenth problem we need on the one hand such time-dependent 
physical (adiabatic) processes. On the other hand, the theory of Quantum Mechanics 
can be used to identify the ground state through the usual statistical predictions from 
the Schrodinger equation with a finitely truncated number of energy states of the time- 
dependent Hamiltonian 7i(t/T). This way, we can overcome the problem of which states 
are to be included in the truncated basis for a numerical study of Quantum Mechanics. 
This also reconciles with the Cantor diagonal arguments which state that the problem 
could not be solved entirely in the framework of classical computation. 

Below is an algorithm [17] based on this philosophy of exploiting the interplay between 
the presumably infinite physical world and the theory of Quantum Mechanics calculated 
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in a finite manner on Turing machines. The algorithm presented may not be the most 
efficient; there could be many other variations making better use of the same philosophy 

It is in general easier to implement some Hamiltonian than to obtain its ground state. 
We thus should start the computation in yet a different and readily obtainable initial ground 
state, \gi), of some initial Hamiltonian, Hj, then deform this Hamiltonian in a time T into 
the Hamiltonian whose ground state is the desired one, through a time-dependent process 
represented by the interpolating Hamiltonian 7i{t/T). 

In this approach, inspired by the quantum adiabatic approach, one starts, for example, 
with a Hamiltonian Hi, 

K 

H i = XX a l ~~ a *i )( a i ~ a 0> ( 15 ) 

which admits the readily achievable coherent state \gj) = \ci\ ■ ■ - ax) as the ground state. 
Then, one forms the time-dependent Hamiltonian 1~L{t/T) in (9), which interpolates in the 
time interval t e [0, T] between the initial Hj and H P . 

• Step 0: Choose an evolution time T, a probability p which can be made arbitrarily 
closed to unity, and an accuracy < e < 1 which can be made arbitrarily small. 

• Step 1 (on the physical apparatus): Perform the physical quantum time- dependent 
process which is governed by the time-dependent Hamiltonian 7i(t/T) and terminates 
after a time T. Then, by projective measurement (either of the observable H p or 
the number operators {N ll . . . , N K }) we obtain some state of the form | • • -rij • • •), 
i = l,--.,K. 

• Step 2 (on the physical apparatus) : Repeat the physical process in Step 1 a number of 
times, L(e,p), to build up a histogram of measurement frequencies (for all the states 
obtained by measurement) until we get a probability distribution P(T; e) at the time 
T with an accuracy e for all the measured states. The convergence of this repetition 
process is ensured by the Weak Law of Large Numbers in probability theory [18]. 
(An overestimate of the number of repetitions is L > l/(e 2 (l — p)).) Note the lowest 
energy state so obtained, \n c ), as the candidate ground state. 

• Step 3 (on the classical computer) : Choose a truncated basis of M vectors made up of 
\ai ■ ■ - ax) and its excited states by successive applications of the displaced creation 
operators b\ = (a\ — a*) on the initial state. 

• Step 4 (on the classical computer): Solve the Schrodinger equation in this basis for 
-0(T), with the initial state ip(0) = \cti • • - ax), to derive a probability distribution 
P cst (T; M) (through \ {^{T)\ ■■■n i ---)\ 2 ) which is similar to that of Step 2 and which 
depends on the total number M of vectors in the truncated basis. 

• Step 5 (on the classical computer): If the two probability distributions are not uni- 
formly within the desired accuracy, that is, \P est (T; M) — P(T;e)| > e, we enlarge 
the truncated basis by increasing the size M and go back to the Step 4 above. 
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Step 6 (on the classical computer): If the two probability distributions are uniformly 
within the desired accuracy, that is, |P est (T; M)—P(T; e)| < e, then use this truncated 
basis to diagonalise Hp to yield, within an accuracy which can be determined from 
e, the approximated ground state \g') and its energy E g >. 

Step 7 (on the classical computer): We can now estimate in this truncated basis the 
gap between the ground state and the first excited state. From this gap, we can make 
use of the quantum adiabatic theorem and choose a time T such that the system has 
a high probability to be in the ground state 



IW(T)>| 2 -1 



< e. 



We then go back to Step 1 with this choice of T, which is to amplify and thus confirm 
the candidate ground state as the real ground state. 

Our point is on computability and not on computational complexity, which depends on 
individual polynomials. Computability is based on the arguments that the adiabatic time 
T is finite [19, 20] (for a high probability of achieving the ground state) and that the 
ground state can be verified by employing the theory of Quantum Mechanics. As long as 
the energy gap is finite so is the computational time. In contrast, the most general classical 
algorithm for Hilbert's tenth problem (by systematic substituting in integers of larger and 
larger magnitudes) cannot solve it in principle even allowing for exponentially grown, but 
finite, amount of time - unless infinite amount of time were available, which it is not. 

Given a Diophantine equation, the substitution in integers of larger and larger magni- 
tudes is not satisfactory as we do not know when the substitution should be terminated. 
Likewise, if we want to numerically simulate the quantum algorithm proposed, we would 
have to use a finitely truncated having M vectors. But we face the same problem of not 
knowing which M to choose in general. That is why the problem is noncomputable on 
classical computers. 

Even we can estimate T, as in the appendix, for some starting point of the adiabatic 
computation as we might not be able in general to exactly know T a priori. Were the re- 
quired adiabatic time T somehow known exactly within classical computation and without 
the help of quantum computation then the problem might be solved classically. But as 
Hilbert's tenth problem cannot be solved by classical computation, we will have to resort to 
quantum computation without a priori knowing exactly the time T, except the knowledge 
that it is finite. Constructive logicians [21] allow for this algorithmic situation under the 
so-called Markov's Principle. 



9 Discussion of the algorithm 

The quantum algorithm above can be proved to terminate (even though it could be after 
a very long time) and give us the decision result for Hilbert's tenth problem. 

The real spectrum of H p is of integer values (in suitable units), and that is what we 
also get from measurement. But the spectrum calculated from a finitely truncated basis 
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is not of integer values and will fluctuate with fluctuation size depends on the size of the 
truncated basis employed. The accuracy size e of the measured probability distribution 
is chosen such that the off-set of the ground state energy 5 should allow us to conclude 
whether the ground state energy E g < is zero or not. (5 is in general a function of e and T.) 

Eg, = {g'\H P \g'), 

= E c + (r\H P \r) + 2Re(r\H P \n c ), 

= E c +(r\H P \r) + 2E c (Re(r\n c }), 

= E c + 5(e), (16) 

where |r) = \g') — \n c ) and Hp\n c ) = E c \n c ). 

The termination of our algorithm is obtained if and when the adding of higher 6-number 
states (those created from the coherent state by the application of the creation operator 
b\ = (aj — a*)) to the truncated basis does not change the approximated ground state of 
Hp beyond certain range of accuracy, 

|<5(e) | < E c ^0. (17) 

Even we can prove that the approximated ground state \g') and its energy will eventually 
converge to its true values, the mathematical noncomputability results from the fact that 
their rates of convergence are unknown. Thus, we might not be able to use mathematical 
reasoning alone to determine when to stop adding more states to the truncated basis in 
order to approximate the ground state correctly. Different truncated bases would give some 
estimates for the ground state but we have no control over these estimates and no idea 
how good they are. They could be anywhere in relation to the true values. This is nothing 
but mathematical noncomputability. (However, this quantum algorithm has inspired us to 
reformulate Hilbert's tenth problem with continuous variables [20], and the mathematical 
computibility of this reformulation, or lack of it, should be investigated further.) 

To know when the truncated basis is sufficiently large to have the estimated ground 
state values within any given accuracy, that is, to regain computability, we have to exploit 
the measurability of physical processes. Because of this measurability we can estimate the 
true accuracy of our measured values. Then a comparison of results from the Schrodinger 
equation to these measurable quantities will help determining the accuracy of results from 
the equation, that is, regaining the lost computability through the physical world, presum- 
ably infinite. 

10 Decoherence and error correction 

Our approach above is in fact a combination of the quantum computation of continu- 
ous variables and of adiabatic evolution. There exists some error correction protocol for 
continuous variables [22] which could be of help here to protect the wave functions from 
decoherence. However, the adiabatic computation we exploit is quite robust in general [23]. 
An imperfect conventional quantum algorithm might have different sorts of errors than an 
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imperfect adiabatic process, where the system is kept close to the instantaneous ground 
state over time. Decoherence by the environment inducing transition between levels could 
be controlled in principle at a temperature that is small compared to the gap (13), given 
its estimate. Errors introduced by the hamiltonian simulation may result in a hamiltonian 
different from (9) but in a form H(t/T) + K(t). Recent numerical study [23] of small 
systems has in fact indicated that the adiabatic computation is interestingly robust even 
for fairly large K(t), provided K(t) varies either sufficiently slowly or sufficiently rapidly 
(which is more likely to be the case considered here due to the nature of our hamiltonian 
simulations). 

11 Probability boosting 

Due to the particular nature of the present scheme, the various approximations result in 
some probability that the final measurement will not find the system in the desired ground 
state; but appropriate choices of the various time parameters could increase the success 
probability of the algorithm to more than even. We spell out explicitly here the probability 
boosting technique, as mentioned in [5], that one could subsequently apply. 

In our computation, adiabatic or otherwise, the end result, starting from some initial 
state \gi), may be contaminated with some excited states other than the desirable ground 
state \g), 

\ 9I ) a\g) + b\e). (18) 

If \a\ > \b\, that is, if there is a better-than-even chance to obtain \g) by measurement 
then one can boost the success probability to arbitrarily closed to unity by performing a 
concatenated computation over I Hilbert spaces 

i 

\gi)i®---\gi)i ^ NY. aPhl ~ p \z)i®---®\g)k 1 ®---®\g)k v ®---®\e) h (19) 

p=0 

where M is the normalising factor. Let us consider the majority amplitudes, when more 
than half of the / Hilbert spaces return the correct results, J\fa p b l ~ p ,p > |; and the minority 
amplitudes, Afa q b l ~ q , q < |. The ratio of the majority over the minority, (a/b) p ~ q , is clearly 
boosted for sufficiently large / and for p ~ 0(1) and q ~ 0(1). The probability distribution 
for a measurement of the end state in (19), as a consequence, is exponentially dominated 
by the majority results. In other words, the probability to obtain a majority result which 
contains the true ground state can be made arbitrarily closed to unity, (1 — e'), provided 
\a\ > \b\ and I > — Cloge'. However, the decoherence control for such I concatenated 
Hilbert spaces will be more crucial. 

12 Concluding remarks 

In this paper, we consider and emphasise on the issue of computability in principle, not 
that of computational complexity. This attempt of broadening of the concept of effective 
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computability, taken into account the quantum mechanical principles, has been argued to 
be able in principle to decide the classically /Turing undecidables, Hilbert's tenth problem 
and thus the Turing halting problem in this instance. If this is realisable, and we don't 
have any evidence of fundamental nature to the contrary, the Church- Turing thesis should 
be modified accordingly. 

In summary, we have encoded the answer to the question about the existence or lack 
of non-negative integer solutions for any Diophantine equation into the of ground state of 
some relevant Hamiltonian. The key factor in the ground state verification is the prob- 
ability distributions, which not only can be calculated in numerical Quantum Mechanics 
(with a truncated basis) but also are measurable in practice. After all, probability dis- 
tributions are also physical observables. However, in using the probability distributions 
as the identification criteria, we have to assume that Quantum Mechanics is able to de- 
scribe Nature correctly to the precision required. Note also that we have here a peculiar 
situation in which the computational complexity, that is, the computation time, might 
not be known exactly before carrying out the quantum computation - although it can be 
estimated approximately (see the appendix below). 

On the other hand, if for any reasons the algorithm is not implementable because of 
physical principles and/or physical resources then it would be an example of information 
being limited by physics, rather than by logical arguments alone. 

Our study is an illustration of "Information is physical" (see [24] for another quantum 
mechanical approach and [25] for where the theory of General Relativity is also exploited 
for the computation of Turing noncomputables). 

That some generalisation of the notion of computation could help solving the previous 
undecidability/noncomputability has been recognised in mathematics and was considered 
by Kleene as quoted in [4]. But this has not been realisable until now simply because of 
the non-recognition of quantum physics as the missing ingredient. Our quantum algorithm 
could in fact be regarded as an infinite search through the integers in a finite amount of 
time, the type of search required by Kleene to solve the Turing halting problem. 

Our decidability study here only deals with the property of being Diophantine, which 
does not cover the property of being arithmetic in general, and as such has no direct 
consequences on the Godel's Incompleteness theorem [26]. However, it is conceivable that 
the Godel's theorem may lose its restrictive power once the concept of proof is suitably 
generalised with quantum principles. 

Note added in proof 

The author has recently obtained a criterion based on the final propbability distribution, 
which can be measured to arbitrary precision, for identifying the ground state of Hp at 
time T. See [28] for more details. 
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Appendix: Gap estimation 

The question of computational complexity, i.e. how large the adiabatic computational time 
T is for a high probability of measurement success, is dependent on Hp, i.e. on the specific 
Diophantine equation in question, and on the initial hamiltonian Hi. Some estimate for 
the energy gap (13) is desirable, but not necessary as discussed above. In this appendix 
we propose such an estimate for the Step of the quantum algorithm. 

It has been shown elsewhere [19, 20] that in general there should be no level crossing for 
the ground state except at the end points t — 0, T where the adiabatic process can start or 
end with some obvious symmetry. Furthermore, the freedom in choosing the initial hamil- 
tonians Hi and their ground states and in performing different adiabatic interpolations, 
not just as linear as in (9), might be exploited to enable the gap enlargement and to speed 
up the computation. 

We shall employ the simple harmonic oscillators, i.e. gaussian approximations, to obtain 
an estimate for the energy gap. 

For the various di, a] appearing in the interpolating hamiltonian Ti{s) (9) at the time 
instant s, we use the Bogoliubov ansatz, with real Uj(s) and Vi(s), 

Ci = Ui(s)ai + Vi(s)a\, 

at = Ui{s)bi - Vi{s)c\. (i) 

The Q-zero occupation state is denoted by |0 C ), where the s-dependence of the state is 
implicit. The canonicity [q,c]] = Sij demands 

«?(*)-«?(*) = 1, (ii) 

upon which 

(a\aj) = vf(s)Sij, 

(a\a)) = (aidj) = -^(s)^(s)<^-, (iii) 

where 

(A k .A k .) = (0cl.4fc.i4fc. 1 C ), (iv) 
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We pay particular attention to the following terms 

: ajdi : c = (u?(s) + v 2 (s)) cJq - Ui(s)vi(s) (cf 2 + c 2 ) , 

: didi : c = -2^(s)q(s)cfq + v 2 (s)cf 2 + m 2 (s)c- . (v) 

Needed next is a version of the Wick theorem [27] for the ordinary product involving the 
operators A 1 , A n , 

Ai - • • A n = ^ : A\ ■ ■ ■ A kl ■ ■ ■ Ak 2p ■ ■ ■ A n : c 

x {(A kl A k2 ) ■ ■ ■ (Ak^A,,^} + permutations} , (vi) 

where the normal ordering : • • • : c is done with respect to some annihilation operators q 
which annihilate some common state |0 C ) for various i, Cj|0{,) = 0. The equation (vi), in 
which the hatted operators are omitted from the normal ordering, is an exact result and 
can be proved by induction. 

We can now list the various steps of our estimation: 

1. We apply the Wick theorem, with respect to q's, to the hamiltonian H(s) (9) 

H(s) = E b (s)I 

+ ^2 terms linear in q and cf 



+£ (Gi(s)cfe + ^( S )c 2 + ^;( s )cf 2 ) 



+ terms involving cfcj, CiCj, cfc] 

+ ^] higher-order, normal-ordered terms of c and c^. (vii) 

ijk— 

Note that in this process the higher-order terms contribute to the coefficients of 
lower-order ones through products of various expectation values in (iii). With help 
from (iii) and (v), the various coefficients Ej,(s), Gi(s), Ki(s), ... in the right hand 
side above can be expressed as polynomials in (wj(s), v i(s)). 

2. We next fix Ui(s), and q(s) numerically from (ii) and from the imposition that the 
coefficients K^s) = in (vii). Then we can evaluate the coefficients Gj(s) from their 
polynomial expressions in (ui(s) , Vi(s)) . From (vii), on the other hand, we see also 
that 

Gi(s) = (l Ci \H(s)\l Ci )-(0 c \H(s)\0 c ), (viii) 

where |l Ci ) = cf|0 c ). Mathematically, minj|Grj(s)| thus provides some indication of 
the size of the energy gaps of Ti(s) around the energy level E b (s) = (0 c \H(s)\0 c ). Even 
though |0 C ), obtained from the linear Bogoliubov transformations (i), in general may 
not be the true ground state of 7~C{s), we could use minj|Gj(s)| as some indicator for 
the gap g in (13) in order to estimate T from the adiabaticity condition (11). 
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Some variations of the above method might yield a better estimate. For instance, one 
could numerically obtain (ui(s), i>j(s)), and thus minj|Gj(s)|, by minimising the energy 
Eb(s) subjected to the constraints (ii) with or without the constraints Ki(s) = 0. The 
aim of this minimisation, if possible, is to select a state |0 C ) whose energy expectation 
value at the time s is as closed to that of the ground state at that instant as allowed by 
the linear Bogoliubov transformations. Other variations might be involving, instead of (i), 
some non-linear relations, for the canonicity only requires that a and c are to be unitarily 
transformed, c = U^(a^,a)aU(a^ ,a). 

The accuracy of the estimation and its higher order corrections can be evaluated sys- 
tematically from the higher order terms in the last line of (vii). Numerical diagonalisation 
of (vii) with a series of truncated bases in \n Ci ) of increasing sizes could further provide us 
more information about the gap thus obtained. 
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